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Relative  Contributions  of  Inelastic 
and  Elastic  Diffuse  Phonon 
Scattering  to  Thermal  Boundary 
Conductance  Across  Solid 
Interfaces 

The  accuracy  of  predictions  of  phonon  thermal  boundary  conductance  using  traditional 
models  such  as  the  diffuse  mismatch  model  (DMM)  varies  depending  on  the  types  of 
material  comprising  the  interface.  The  DMM  assumes  that  phonons,  undergoing  diffuse 
scattering  events,  are  elastically  scattered,  which  drives  the  energy  conductance  across 
the  interface.  It  has  been  shown  that  at  relatively  high  temperatures  (i.e.,  above  the 
Debye  temperature)  previously  ignored  inelastic  scattering  events  can  contribute  sub¬ 
stantially  to  interfacial  transport.  In  this  case,  the  predictions  from  the  DMM  become 
highly  inaccurate.  In  this  paper,  the  effects  of  inelastic  scattering  on  thermal  boundary 
conductance  at  metal/dielectric  interfaces  are  studied.  Experimental  transient  thermore¬ 
flectance  data  showing  inelastic  trends  are  reviewed  and  compared  to  traditional  models. 
Using  the  physical  assumptions  in  the  traditional  models  and  experimental  data,  the 
relative  contributions  of  inelastic  and  elastic  scattering  to  thermal  boundary  conductance 
are  inferred.  [DOI:  10.1115/1.2995623] 

Keywords:  thermal  boundary  conductance,  diffuse  phonon  scattering,  elastic  scattering, 
inelastic  scattering,  nanoscale  solid  interfaces 


Introduction 

An  understanding  of  the  basic  energy  transport  mechanisms 
involved  in  interfacial  thermal  transport  is  critical  for  thermal 
management  of  nanostructured  devices.  When  the  length  scale  of 
a  device  is  comparable  to  or  shorter  than  the  host  material’s  ther¬ 
mal  diffusion  length,  which  is  often  the  case  for  modern  nanoscale 
devices,  heat  transport  away  from  the  active  regions  is  greatly 
affected  by  the  interfacial  properties  [1],  An  ever  increasing  chal¬ 
lenge  in  the  development  of  these  devices  is  successfully  engi¬ 
neering  the  heat  transport  across  the  interfaces  to  control  the  ther¬ 
mal  responses  experienced  in  the  materials  in  the  devices.  This 
involves  a  fundamental  understanding  of  the  thermal  boundary 
conductance,  hBD,  at  the  specific  interfaces. 

This  study  examines  the  relative  contributions  of  inelastic  and 
elastic  diffuse  phonon  scattering  on  thermal  boundary  conduc¬ 
tance.  Several  models  to  predict  thermal  boundary  conductance 
are  discussed  in  detail  and  compared  to  experimental  data.  Based 
on  the  assumptions  of  these  models,  a  new  model  to  predict  the 
maximum  hBD  due  to  inelastic  scattering  is  presented.  Using  this 
new  model,  the  relative  contributions  of  elastic  and  inelastic  scat¬ 
tering  on  hBD  are  examined. 

Models  for  Phonon  Thermal  Boundary  Conductance 

In  an  attempt  to  predict  thermal  boundary  conductance  at  low 
temperatures.  Little  [2]  proposed  the  acoustic  mismatch  model 
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(AMM)  to  account  for  the  specular  scattering  of  phonons  at  an 
interface  between  two  materials.  In  the  following  derivation  of  the 
AMM  and  subsequent  discussion,  side  1  will  refer  to  the  side  of 
the  interface  with  the  “softer”  material  and  side  2  will  refer  to  the 
side  of  the  interface  with  the  “stiffer”  material.  The  softer  mate¬ 
rials  are  characterized  by  their  lower  phonon  velocities,  smaller 
phonon  vibrational  spectrum,  and  lower  Debye  temperatures  com¬ 
pared  to  the  stiffer  materials,  which  have  higher  phonon  veloci¬ 
ties,  larger  phonon  vibrational  spectrum,  and  higher  Debye  tem¬ 
peratures.  For  example,  consider  a  metal/dielectric  system  in 
which  phonons  are  propagating  in  a  metallic  film  such  as  Pb  (the 
lower  Debye  temperature,  softer,  material)  toward  the  interface 
with  a  dielectric  substrate  such  as  diamond  (the  higher  Debye 
temperature,  stiffer,  material).  In  this  case,  the  metallic  Pb  film  is 
referred  to  as  side  1  and  the  dielectric  diamond  substrate  as  side  2. 
Table  1  gives  phonon  velocities,  vibrational  cutoff  frequencies, 
and  Debye  temperatures  for  several  materials  of  interest  in  this 
work.  The  net  heat  flux  occurring  from  side  1  to  side  2  can  be 
calculated  by 

4  =  (1) 

Using  the  analogy  between  photons  and  phonons  as  wavepackets 
of  energy,  the  intensity  of  phonons  [3]  is  used  to  calculate  the  heat 
flux  determined  from  the  equation  of  phonon  radiative  transfer 
(EPRT)  [4],  The  heat  flux  is  expressed  as 

1  (n'2 

4  =  “2  Dlj{(n)n{ix),T)fKvvljalj{(l>,(n) 

1  j  Jo  Jo 

Xcos(<f)sin(<fi)d(od<fi 

=  hBDAT  (2) 

where  coc  is  the  cutoff  frequency,  D(to)  is  the  density  of  states, 
n(co,T )  is  the  Bose-Einstein  phonon  distribution  function,  to  is  the 
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Table  1  Debye  temperature  [8,24,25],  longitudinal  and  transverse  phonon  velocities  [19,25], 
and  calculated  cutoff  frequencies  of  materials  of  interest  in  this  study 


Material 

MK) 

vL  (m  s  l) 

vT  (m  s  *) 

ofL  X  10"13  (s-1) 

u>cTX  10-13  (s-1) 

A1203 

1035 

10,890 

6450 

12.1 

7.19 

AIN 

1150 

11,120 

6267 

15.8 

8.89 

Au 

165 

3390 

1290 

5.14 

1.96 

Bi 

119 

1543 

1107 

1.45 

1.04 

Diamond 

2230 

17,500 

12,800 

30.2 

22.1 

Pb 

105 

2350 

970 

2.96 

1.22 

Pt 

240 

4174 

1750 

6.60 

2.77 

phonon  frequency,  v  is  the  phonon  velocity,  and  a((f>,co)  is  the 
phonon  transmission  probability.  In  this  analysis,  a  local  thermo¬ 
dynamic  equilibrium  is  assumed  around  the  phonon  scattering 
event,  so  local  temperatures  can  be  defined.  The  subscripts  “1” 
and  “/’  refer  to  the  side  and  the  phonon  mode  (longitudinal  or 
transverse),  respectively.  Rearranging  Eq.  (2)  and  for  the  limit  of 
AT  approaching  zero,  the  AMM  is  given  by 


-  j  Jo  Jo 


'w  dn((o,  T) 

Di  j(") — — — fiojvtjotfjc/),  id) 


Xcos((f))sm((l>)dwd(f>  (3) 

Note,  in  the  AMM,  the  transmission  probability  is  a  function  of 
incident  angle.  Although  Eq.  (3)  has  been  shown  to  predict  hBD 
relatively  well  at  low  temperatures  (T<  7  K)  and  at  ideal  inter¬ 
faces  where  specular  scattering  is  probable  [5],  this  represents 
only  a  very  limited  population  of  interfaces  in  modern  devices, 
which  may  operate  at  higher  temperatures  and  have  disordered 
regions  near  the  interface  that  induce  diffuse  scattering. 

To  account  for  this  type  of  phonon  scattering,  Swartz  and  Pohl 
[6]  developed  the  diffuse  mismatch  model  (DMM)  to  predict  hBD 
at  more  realistic  interfaces.  The  DMM  theory  assumes  elastic 
scattering,  that  is,  upon  scattering  a  phonon  forgets  where  it  came 
from  and  thus  the  probability  of  reflection  from  one  side  equals 
the  probability  of  transmission  from  the  other.  The  probability  that 
a  phonon  will  scatter  into  a  given  side  of  the  interface  is  thus 
independent  of  where  it  came  from,  but  rather  it  is  proportional  to 
the  density  of  phonon  states  on  that  side  and  is  restricted  by  the 
principle  of  detailed  balance. 

To  apply  the  DMM  in  its  simplest  form,  the  following  assump¬ 
tions  must  be  made  [7]:  (1)  phonons  are  elastically  scattered,  i.e., 
a  phonon  from  side  1  with  frequency  co  can  only  emit  a  phonon 
from  the  interface  with  the  same  frequency  co;  (2)  phonon  scatter¬ 
ing  is  completely  diffuse,  i.e.,  a  scattered  phonon  has  no  memory 
of  the  mode  (longitudinal  or  transverse)  or  direction  of  the  inci¬ 
dent  phonon;  and  (3)  the  materials  on  both  sides  of  the  interface 
are  elastically  isotropic,  i.e.,  the  longitudinal  and  transverse  acous¬ 
tic  velocities  are  constant  in  all  crystallographic  directions.  As¬ 
sumption  (3)  relaxes  the  angle  dependence  in  Eq.  (3),  which  leads 
to 

hBDM=  ^2  Dj  f  1  a^hcoDijico)  dw  (4) 

j  Jo 

The  linear  Debye  approximation  for  the  phonon  dispersion  will  be 
used  to  calculate  the  phonon  density  of  states.  The  cutoff  fre¬ 
quency  of  each  phonon  mode  j  on  side  1  can  be  calculated  for  the 
separate  phonon  modes.  Assuming  Debye  dispersion,  the  cutoff 
frequency  can  be  calculated  by  coj  ■=  f1>/-(6'7T2A’i)1/3,  where  N\  is 
the  total  number  of  oscillators  per  until  volume  of  side  1  [8].  In 
cubic  structures  (such  as  metal  with  one  atom  per  unit  cell),  N\  is 
simply  the  atomic  density,  calculated  by  Ni=pNA/M,  where  p  is 
the  mass  density,  NA  is  Avogadro’s  number,  and  M  is  the  atomic 
weight.  However,  in  structures  with  more  than  one  atom  per  unit 
cell  (for  example,  diamond  structures  with  diatomic  basis  such  as 


Si  or  diamond),  the  number  of  primitive  cells  per  unit  volume 
must  be  divided  by  the  number  of  atoms  in  the  basis  [9].  There¬ 
fore,  for  a  diatomic  basis,  Nl=pNAl{2M). 

Due  to  the  assumption  of  diffuse  scattering,  the  probability  of 
transmission  from  side  1  to  side  2  is  the  same  as  the  probability  of 
reflection  from  side  2  to  side  1,  i.e.,  ai{co)  =  \- a2(to).  Therefore, 
according  to  the  principle  of  detailed  balance,  the  phonon  trans¬ 
mission  probability  is  calculated  by  [6] 


(5) 


This  simplified  transmission  coefficient  is  a  result  of  assuming  a 
Debye  density  of  states  and  elastic  scattering.  The  transmission 
probability  calculated  with  the  Debye  density  of  states  agrees  well 
with  the  transmission  of  two  simple  structures  calculated  with  a 
more  realistic  density  of  states  and  molecular  dynamics  simula¬ 
tions  (MDSs)  [10,11],  The  DMM  has  been  shown  to  predict  the 
response  of  higher  temperature  interfaces  (T>  15  K)  relatively 
well  [5,6].  However,  at  much  higher  temperatures,  the  DMM  has 
been  shown  to  do  a  poor  job  of  predicting  the  value  of  hBD,  in 
some  cases  underpredicting  [12]  while  in  other  cases  overpredict¬ 
ing  [7,12-14], 

Another  model  for  hBD,  the  phonon  radiation  limit  (PRL),  esti¬ 
mates  the  maximum  conductance  for  interfacial  transport  due  to 
elastic  scattering  [15].  The  development  of  the  PRL  is  very  similar 
to  that  of  the  DMM.  However,  the  PRL  assumes  that  all  of  the 
phonons  in  side  2  up  below  the  cutoff  frequency  in  side  1  (assum¬ 
ing  &)[  <  co2 )  contribute  to  thermal  transport  though  elastic  colli¬ 
sions  (o'=l).  With  these  assumptions,  the  PRL  is  given  as 


h 


PRL  _ 
BD  ~ 


hcoDifco) 


dn(co,T ) 
dT 


dco 


(6) 


which  represents  the  maximum  conductance  due  to  elastic  scatter¬ 
ing. 

The  applicability  of  all  of  the  aforementioned  models  is  limited 
to  interfaces  in  which  thermal  transport  is  dominated  by  elastic 
phonon  scattering.  Inelastic  phonon  scattering  has  been  shown  to 
offer  an  additional  channel  for  thermal  transport,  which  can  lead 
to  different  values  and  trends  than  predicted  with  available  models 
[10,16,17],  In  the  classical  limit  (for  real  materials  T>  0d,u  the 
Debye  temperature  of  side  1,  the  lower  Debye  temperature  mate¬ 
rial),  hBD  calculated  by  either  the  DMM  or  PRL  is  relatively  in¬ 
dependent  of  temperature.  The  only  temperature  dependent  part  of 
either  model  is  in  the  distribution  function,  which  at  temperatures 
well  above  the  Debye  temperature  becomes  constant.  Because 
these  models  do  not  assume  any  inelastic  scattering,  the  models 
are  independent  of  temperature  at  high  temperatures  (when  T 
>  0D  i).  To  check  the  temperature  dependence  of  hBD,  Stevens  et 
al.  [10]  conducted  several  molecular  dynamic  simulations  at  a 
range  of  temperatures.  A  strong  linear  relationship  was  observed 
in  the  results  of  the  MD  calculations.  This  linear  trend  in  hBD 
when  T>  ffDl  has  been  observed  experimentally  with  transient 
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Fig.  1  DMM  (solid  line),  PRL  (dotted  line),  and  JFDMM  (dashed 
line)  calculations  compared  to  temperature  dependent  TTR 
data  on  (a)  Pt/AIN  [17],  ( b )  Pt/Al203  [17],  (c)  Au/diamond  [14], 
( d)  Bi/H/diamond  [16],  and  (e)  Pb/diamond  and  Pb/H/diamond 
[16].  Note  the  similar  temperature  dependent  trends  of  the 
DMM  and  PRL  due  to  their  assumption  of  elastic  scattering. 


thermoreflectance  (TTR)  pump-probe  measurements.  Stoner  and 
Maris  [14]  have  shown  this  linear  trend  for  a  Au  film  on  a  dia¬ 
mond  substrate.  This  trend  was  shown  on  several  other  Him  sub¬ 
strate  material  systems  at  low  temperatures  by  Lyeo  and  Cahill 
[16]  and  at  high  temperatures  by  Hopkins  et  al.  [17],  In  all  cases, 
values  and  trends  in  the  experimental  data  differ  from  those  pre¬ 
dicted  by  the  DMM  and  PRL. 

An  underlying  assumption  governing  the  DMM  and  PRL  is  that 
a  phonon  transmits  energy  across  an  interface  by  emitting  a  pho¬ 
non  with  the  same  frequency,  i.e.,  the  phonons  are  elastically  scat¬ 
tered.  Therefore,  as  the  interface  temperature  is  driven  up  above 
the  Debye  temperature  of  the  softer  material,  hBD  is  predicted  to 
be  relatively  constant  by  the  DMM.  For  example,  Pb/diamond  in 
Fig.  1  where  0/3Pb=lO5  K  and  0D  diamond=2200  K  [8].  This  is  a 
result  of  the  change  in  phonon  population  with  temperature  pre¬ 
dicted  by  the  Bose-Einstein  distribution  function.  At  temperatures 
close  to  a  material’s  Debye  temperature,  the  change  in  phonon 
population  with  temperature  becomes  linear.  The  DMM  is  a  func¬ 
tion  of  the  temperature  derivative  of  the  Bose-Einstein  distribu¬ 
tion,  which  results  in  the  constant  hBD  predicted  at  higher  tem¬ 
peratures  (T>  6d  Pb).  Assuming  elastic  scattering,  the  hBD 
predictions  would  follow  a  trend  relating  to  the  change  in  the  Pb 
phonon  population  with  temperature.  However,  if  inelastic  phonon 
processes  occur  (i.e.,  a  phonon  with  frequency  scatters 

into  several  phonons  with  frequencies  below  &)pb),  then  the  change 
in  hBD  with  temperature  would  be  related  to  the  change  in  the 
diamond  phonon  population  in  addition  to  the  Pb  population. 

Using  this  approach,  Hopkins  and  Norris  [18]  developed  a 
simple  correction  to  the  DMM  to  account  for  the  discrepancy 
between  the  DMM  and  the  experimental  data  in  the  event  of  in¬ 


elastic  scattering.  By  blending  the  vibrational  spectra  of  the  film 
and  substrate  materials,  an  approximation  was  developed  for  the 
contribution  of  these  inelastic  modes  with  a  diffuse  scattering  as¬ 
sumption,  the  joint  frequency  diffuse  mismatch  model  (JFDMM). 
The  JFDMM  assumes  the  same  form  as  the  DMM  (Eq.  (4))  but 
uses  a  modified  phonon  velocity  that  is  taken  as  a  weighted  aver¬ 
age  of  the  velocities  of  the  phonons  of  sides  1  and  2.  Conse¬ 
quently,  this  results  in  a  weighted  average  of  the  phonon  spectra 
used  in  the  hBD  calculation,  given  by 

Dmod,j  ~  i  3  ’  ra  ^  wmodj  (^) 

—  ^  Lnodj 


=  W61 rt&Ni  +  &V2)}1/3 

(8) 

*,modj=fl*T  +  &»'2 

(9) 

where  the  weighting  factor  f  is  simply  a  percentage  of  the  com¬ 
position  of  each  material  in  the  unit  volume,  mathematically  ex¬ 
pressed  as 


£t  = 


IVi 

—M, 

N2 


—M{  +  M 

N2 


1 


(10) 


where  M  is  the  atomic  mass.  This  approximation  introduces  high 
frequency  phonons  that  are  available  in  the  vibrational  spectrum 
in  side  2  but  not  side  1  into  the  incident  heat  flux.  The  JFDMM 
increases  the  prediction  of  hBD  by  a  factor  that  is  proportional  to 
the  side  2  vibrational  spectrum,  giving  an  approximation  for  in¬ 
elastic  scattering.  The  transmission  coefficient  for  the  JFDMM  is 
still  assumed  to  be  calculated  with  Eq.  (5);  note  that  this  relaxes 
the  assumption  of  detailed  balance  since  the  JFDMM  assumes  a 
modified  phonon  flux. 

Figure  1  compares  the  predictions  from  the  DMM,  PRL,  and 
JFDMM  to  experimental  TTR  data  taken  on  several  different  ma¬ 
terial  systems  over  a  wide  range  of  temperatures.  This  figure  com¬ 
pares  predictive  trends  to  data  on  (a)  Pt/AIN  (Hopkins  et  al.)  [17], 
(b)  Pt/Al203  (Hopkins  et  al.)  [17],  (c)  Au/diamond  (Stoner  and 
Maris)  [14],  (d)  Bi/H/diamond  (Lyeo  and  Cahill)  [16],  and  (e) 
Pb/diamond  and  Pb/H/diamond  (Lyeo  and  Cahill)  [16],  The 
model  calculations  use  elastic  constants  to  calculate  phonon  ve¬ 
locities  and  material  properties  to  calculate  the  cutoff  frequencies 
[19],  Note  that  the  models  cannot  distinguish  between  different 
deposition  or  interface  conditions  [20],  so  the  predictions  by  the 
models  on  the  hydrogen  terminated  and  non-hydrogen-terminated 
substrates  are  the  same.  In  Fig.  1,  the  DMM  is  represented  by  the 
solid  line,  the  PRL  is  represented  by  the  dotted  line,  and  the 
JFDMM  is  represented  by  the  dashed  line.  Note  that  the  JFDMM 
predicts  a  closer  value  and  better  temperature  dependent  trend  to 
the  experimental  data  than  the  DMM,  and  in  most  cases  the  PRL. 
Since  the  DMM  and  PRL  both  assume  elastic  scattering  of  fre¬ 
quencies  only  up  to  the  side  1  cutoff  frequency,  these  models 
share  the  same  trend  with  temperature.  The  JFDMM,  however, 
shows  a  different  temperature  dependent  trend  that  is  more  in  line 
with  the  experimental  data  since  it  assumes  phonons  with  frequen¬ 
cies  higher  than  the  cutoff  frequency  in  side  1  can  participate  in 
hBD.  These  data  represent  material  systems  that  show  some  evi¬ 
dence  of  inelastic  scattering  at  these  temperatures. 


Inelastic  Phonon  Radiation  Limit 

Although  the  JFDMM  shows  improvement  in  hBD  predictions 
in  the  event  that  inelastic  phonon  scattering  dominates  interfacial 
transport,  this  method  makes  necessary  assumptions  about  phonon 
transport  that  deserve  further  attention.  The  main  assumption  of 
the  JFDMM  is  that  a  fraction  of  the  available  phonon  states  in  the 
substrate  are  present  in  the  film.  In  actuality,  the  atoms  around  the 
interface  are  vibrating  at  joint  modes  by  Newton’s  law  of  motion 
[21],  The  maximum  allowable  frequency  for  these  joint  modes  is 
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the  substrate  cutoff  frequency.  However,  dampening  of  these 
modes  will  occur  due  to  the  differing  interatomic  forces  of  the 
film  and  substrate.  This  prevents  the  joint  modes  from  oscillating 
at  certain  frequencies.  This  continuum  phenomenon  is  paralleled 
quantum  mechanically  in  the  JFDMM  through  the  weighted  aver¬ 
age  of  two  materials’  phonon  states  in  the  incident  flux.  The 
weighing  factor,  f,  reflects  the  dampening  of  the  substrate  modes, 
and  is  used  to  enhance  the  incident  phonon  flux  to  take  into  ac¬ 
count  the  joint  vibrations  that  are  allowed  after  dampening  effects. 
However,  without  computationally  expensive  computer  simula¬ 
tions  or  a  rigorous  theoretical  treatment,  this  weighting  factor  can¬ 
not  be  explicitly  determined  for  every  phonon  mode,  and  therefore 
should  be  viewed  as  an  estimation  of  the  joint  modes  participating 
in  hBD.  That  being  said,  the  JFDMM  can  be  viewed  as  a  starting 
point  for  estimating  the  maximum  hBD  due  to  inelastic  scattering. 

Consider  atoms  vibrating  at  joint  modes  around  the  interface 
with  no  dampening  effects.  Therefore,  these  atoms  can  vibrate  at 
all  allowable  frequencies  up  to  the  maximum  allowed  frequency 
in  the  substrate.  Atoms  on  side  1  and  side  2  will  be  coupled  in 
joint  vibrational  modes  with  frequencies  up  to  the  cutoff  fre¬ 
quency  on  side  2.  This  is  paralleled  in  the  quantum  treatment  by 
forcing  fi=0  and  £2=1,  which  conceptually  is  treating  hBD  as  a 
function  of  the  incident  phonon  flux  from  side  2  transmitted  into 
side  1.  This  simplifies  Eqs.  (7)— (9)  to  Dmoij=D2j,  &>mod,/=fc,2j> 
and  vmoi  j  =  v2 j,  respectively,  and  redefines  Eq.  (4)  in  terms  of  the 
flux  transmitted  from  side  2  to  side  1,  given  as 

Kd  =  ^2  *2j  f  2J  ^A(T)ho)D2J(o))^^do)  (11) 

j  Jo 

From  the  nature  of  diffuse  scattering,  the  inelastic  transmission 
probability  is  o?2A(T)  =  1  -a,l1nel(r),  which  is  different  from  the 
elastic  transmission  probability  calculated  with  Eq.  (5).  Equation 
(11)  allows  for  the  possibility  of  higher  frequency  phonons  that  do 
not  exist  in  the  film  to  participate  in  hBD.  Without  knowledge  of 
the  explicit  temperature  dependence  of  a2,  hBD  cannot  be  calcu¬ 
lated.  However,  examining  Eq.  (11)  along  with  experimental  data 
can  give  important  understanding  of  the  role  of  inelastic  phonon 
scattering  in  thermal  boundary  conductance. 

Consider  the  case  where  all  available  substrate  phonons  are 
participating  in  hBD.  In  this  case,  the  probability  that  a  phonon  on 
side  2  is  inelastically  transmitted  (i.e.,  breaks  down  into  lower 
frequency  phonons  and  transmitted  into  side  1)  is  1.  By  letting 
ot2—  1,  Eq.  (11)  becomes  an  expression  for  the  largest  allowable 
thermal  boundary  conductance  due  to  inelastic  scattering,  or  an 
inelastic  phonon  radiation  limit  (IPRL),  expressed  as 
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The  IPRL  assumes  that  all  side  2  phonons  are  transmitted  into  side 
1,  and  does  not  explicitly  take  into  account  elastic  or  inelastic 
scattering  processes.  However,  by  allowing  all  frequencies  of 
phonons  in  side  2  to  transmit  energy  into  side  1,  which  has  the 
lower  cutoff  frequency,  inelastic  scattering  is  implied.  Note  that  in 
this  limit,  similar  to  the  JFDMM,  which  also  takes  into  account 
some  inelastic  scattering,  the  assumption  of  equilibrium  is  relaxed 
and  therefore  the  principle  of  detailed  balance  is  not  enforced. 
Also,  in  the  case  of  an  interface  in  a  homogeneous  material  (i.e., 
an  “imaginary”  interface  between  two  of  the  same  materials),  the 
IPRL  does  not  impose  any  unphysical  interface  resistance  and 
gives  the  same  thermal  flux  across  an  imaginary  interface  in  a 
homogeneous  material  as  that  predicted  by  the  Fourier  law.  The 
other  models  discussed  thus  far  require  accurate  knowledge  of 
interface  transmission  probability  and  the  use  of  a  diffusion- 
transmission  interface  correction  [22]  to  relax  to  flux  predicted  by 
the  Fourier  law  ensuring  thermal  flux  continuity. 

Equation  (12)  is  solely  dependent  on  the  acoustic  properties  of 
side  2,  so  in  the  IPRL,  hBD  on  any  film/substrate  system  only 
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Fig.  2  (Top)  Experimental  measurements  of  hBD  for  various 
interfaces  (see  Fig.  1)  compared  to  their  corresponding  PRL. 
The  PRL  predicts  a  constant  hBD  at  temperatures  above  the 
Debye  temperature  of  the  lower  Debye  temperature  material. 
(Bottom)  Inelastic  phonon  radiation  limit  for  the  four  material 
systems  in  the  top  graph.  The  IPRL  shows  a  linear  increase 
over  a  temperature  range  of  100-400  K,  the  same  trend  that  is 
shown  in  the  experimental  data.  Note,  however,  that  the  slope 
of  the  linear  increase  in  the  IPRL  is  much  greater  than  the  slope 
of  the  linear  increase  in  the  data. 


depends  on  the  stiffer  material.  This  is  apparent  in  Figs.  2  and  3, 
which  compare  the  temperature  dependent  hBD  data  from  Fig.  1  to 
their  respective  IPRL  and  PRL  calculations.  The  experimental 
data  are  graphed  on  the  same  plot  as  their  PRL  calculations.  The 
IPRL  calculations  are  shown  in  the  lower  plots  of  the  figures  since 
the  values  of  the  IPRL  are  orders  of  magnitude  greater  than  the 
data  and  the  PRL.  However,  the  separate  plots  allow  for  easy 
comparisons  of  temperature  trends  between  the  models  and  the 
data. 

The  IPRL  calculations  for  the  three  material  systems  in  Fig.  2 
are  all  identical,  since  the  IPRL  is  only  dependent  on  the  higher 


Fig.  3  (Top)  Experimental  measurements  of  hBD  for  Pt/Al203 
and  Pt/AIN  (see  Fig.  1)  compared  to  their  corresponding  PRL. 
(Bottom)  Inelastic  phonon  radiation  limit  for  the  two  material 
systems  in  the  top  graph.  The  increase  in  hBD  over  the  tempera¬ 
ture  range  of  interest  is  greater  in  the  IPRL  than  in  the  slope  of 
the  linear  increase  in  the  data,  which  is  greater  than  the  in¬ 
crease  in  the  PRL. 
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Fig.  4  Regime  map  of  thermal  boundary  conductance  models 
that  takes  into  account  various  degrees  of  elastic  and  inelastic 
diffuse  phonon  scattering.  The  DMM  and  PRL,  which  take  into 
account  varying  degrees  of  elastic  scattering,  are  paralleled 
with  the  JFDMM  and  IPRL,  which  take  into  account  varying  de¬ 
grees  of  inelastic  scattering. 


the  linear  trend  in  hBD.  Therefore,  the  total  thermal  boundary  con¬ 
ductance  in  the  classical  limit  (T>  Od)  will  be  a  blend  of  both  the 
constant  elastic  and  temperature  dependent  inelastic  phonon  scat¬ 
tering  contributions,  which  can  be  mathematically  expressed  as 

hBD(T)  =  htD  +  K^(T)  (13) 

where  the  superscripts  el  and  inel  denote  the  elastic  and  inelastic 
contributions  to  hBD.  This  is  similar  to  the  phenomenological  ob¬ 
servation  by  Stevens  et  al.  [10]  based  on  MD  simulations.  Equa¬ 
tion  (13)  separates  the  elastic  and  inelastic  contributions  to  hBD  by 
assuming  that  these  two  energy  transfer  mechanisms  can  be 
treated  as  two  thermal  pathways  in  parallel.  Assuming  that  the 
elastic  and  inelastic  scattering  contributions  to  hBD  can  be  sepa¬ 
rated  is  valid  since  in  acoustically  mismatched  materials  such  as 
those  of  interest  in  this  work,  the  number  of  phonons  at  any  avail¬ 
able  frequency,  co,  in  side  1  is  much  greater  than  the  number  of 
phonons  at  that  same  frequency  co  in  side  2,  so  phonons  of  fre¬ 
quencies  co  in  the  side  1  vibrational  spectrum  can  participate  in 
both  elastic  and  inelastic  scattering  events. 

Since  the  PRL  and  the  IPRL  represent  the  upper  limit  to  elastic 
and  inelastic  scattering,  it  is  expected  that  the  elastic  and  inelastic 
contributions  to  hBD  will  be  some  fraction  of  the  PRL  and  IPRL, 
respectively.  Therefore,  Eq.  (13)  can  be  written  more  explicitly  as 
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Debye  temperature  material  (side  2).  The  PRL,  on  the  other  hand, 
depends  on  the  cutoff  frequency  of  the  lower  Debye  temperature 
material  (side  1),  so  therefore  is  specific  for  each  film/substrate 
material  system.  In  all  cases,  the  IPRL  is  greater  than  the  PRL, 
which  is  intuitive  since  the  PRL  assumes  that  only  a  certain  frac¬ 
tion  (up  to  the  film  cutoff  frequency)  of  substrate  phonons  are 
participating  in  hBD  where  the  IPRL  assumes  that  all  available 
substrate  phonons  can  participate  in  hBD.  In  addition,  notice  that 
the  PRL  predicts  hBD  as  a  constant  value  above  the  film's  Debye 
temperature  (similar  to  the  DMM),  where  the  IPRL  predicts  a 
linear  change  of  hBD  with  temperature,  similar  to  the  trends  in  the 
experimental  data.  In  all  cases,  the  linear  change  in  the  IPRL  is 
much  greater  than  the  linear  change  in  the  data.  Figure  3  shows 
the  same  calculations  as  Fig.  2  for  the  Pt/Al203  and  Pt/AIN  data. 
Note  that  the  IRPL  calculations  for  AIN  and  A1203  are  vastly 
different  (bottom),  yet  the  experimental  data  and  PRL  calculations 
for  the  two  samples  are  nearly  identical  (top).  The  measured  De¬ 
bye  temperatures  for  AIN  and  A1203  are  nearly  identical,  which 
would  mean  they  have  similar  maximum  cutoff  frequencies.  How¬ 
ever,  with  a  Debye  approximation,  the  calculations  of  the  cutoff 
frequencies  do  not  match.  This  could  be  a  result  of  assuming  that 
these  solids  can  be  treated  as  Debye  solids,  which  may  not  be 
completely  accurate  for  complex  dielectrics  such  as  AIN  or  A1203. 
However,  for  the  purpose  of  this  analysis,  the  Debye  solid  is  as¬ 
sumed  for  all  materials  since  it  allows  for  a  simple  calculation  to 
elucidate  important  physics  in  interfacial  phonon  scattering  pro¬ 
cesses. 

Figure  4  shows  a  regime  map  of  the  four  different  models  for 
thermal  boundary  conductance  driven  by  diffuse  phonon  scatter¬ 
ing.  The  DMM,  which  takes  into  account  partial  phonon  transmis¬ 
sivity,  is  the  elastic  version  of  the  JFDMM,  which  takes  into  ac¬ 
count  partial  phonon  transmissivity  during  inelastic  scattering. 
The  upper  limits  to  the  DMM  and  JFDMM  models  are  the  PRL 
and  the  IPRL,  respectively,  which  assume  a  phonon  transmissivity 
of  unity  during  elastic  and  inelastic  interfacial  phonon  scattering, 
respectively. 

Relative  Contributions  of  Inelastic  and  Elastic  Phonon 
Scattering 

From  the  trends  in  the  IPRL  as  compared  to  the  PRL  and  ex¬ 
perimental  data,  the  relative  contributions  of  elastic  and  inelastic 
scattering  can  be  examined.  Above  the  film’s  Debye  temperature, 
it  is  apparent  that  contributions  from  elastic  scattering  will  result 
in  a  constant  hBD.  However,  inelastic  scattering  events  will  drive 


hBD(T)  =Ah\ +  B(T)IiBd\T)  (14) 

Where  is  the  PRL  and  /j^l(T)  is  calculated  by  Eq.  (12).  The 
coefficients  A  and  B(T)  are  coefficients  representing  the  fraction 
of  the  maximum  possible  conductance  due  to  each  scattering  pro¬ 
cess  and  will  be  determined  from  the  experimental  data  in  the  next 
section.  Although  A  can  be  estimated  using  Eq.  (5)  so  that  A=1 
-al  =  a2  and  hBD=AhB^=[l -al\hp^'  since  only  elastic  scatter¬ 
ing  is  dealt  with  in  this  case,  the  fundamental  assumptions  driving 
the  evaluation  of  a  are  flawed,  even  in  the  elastic  limit  [23].  For 
example,  in  calculating  Eq.  (5),  equilibrium  is  assumed  (principle 
of  detailed  balance)  although  thermal  transport  is  inherently  a 
nonequilibrium  process.  Also,  in  the  limit  that  both  materials  ad¬ 
jacent  to  the  interface  are  the  same,  Eq.  (5)  becomes  50%.  How¬ 
ever,  in  this  case,  since  there  is  no  difference  in  the  acoustic  prop¬ 
erties,  phonon  transmission  should  be  100%. 

The  temperature  dependency  of  the  inelastic  phonon  transmis¬ 
sion  coefficient,  B(T),  arises  from  the  fact  that  as  temperature 
increases,  there  are  proportionately  more  substrate  phonons  avail¬ 
able  to  break  down  and  scatter  with  lower  frequency  side  1 
phonons.  This  temperature  dependency  is  apparent  from  the  dif¬ 
fering  linear  slopes  of  the  experimental  data  and  the  IPRL  over  the 
temperature  range  of  interest.  Therefore,  Eq.  (11)  can  be  rewritten 
as 

Kd  =  7  2  v2j  f  ’  a2A(T)hcoD2J{co)d  | ’  \co  =  B(T)hB*L 

j  J  0 

(15) 

where  a'2nel(r)=S(T).  It  is  apparent  in  the  data  in  Figs.  2  and  3 
that  the  temperature  dependency  will  be  some  function  of  the 
temperature  dependencies  of  the  phonon  populations  of  both  side 
1  and  side  2.  The  JFDMM  takes  this  into  account  by  considering 
a  weighted  average  of  the  phonon  populations  of  the  two  sides. 
However,  this  correction  results  in  a  temperature  independent  con¬ 
stant  that  enhances  ax,  and  a  temperature  dependent  transmission 
coefficient  should  be  considered  to  understand  the  relative  effects 
of  inelastic  scattering.  Through  the  principle  of  detailed  balance 
invoked  on  the  incoming  phonon  fluxes  from  sides  1  and  2, 
Dames  and  Chen  [9]  developed  a  temperature  dependent  transmis¬ 
sion  coefficient  based  on  the  changes  of  the  phonon  population  of 
the  two  sides,  given  by 
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Fig.  5  Measured  hBD  on  Pb/diamond  with  best  fit  line  extrapo¬ 
lated  to  determine  the  y-intercept  at  3.26  MW  m~2  K_1,  which  is 
assumed  as  the  elastic  contribution  to  hBD  in  the  classical  limit. 
This  value  is  less  than  the  prediction  of  the  PRL,  as  expected, 
but  agrees  well  with  the  predictions  from  the  DMM,  validating 
the  transmission  coefficient  calculations  in  the  DMM  for  a  Pb / 
diamond  interface. 


inel/rrA  J 

“i  =  - 

ZjnjUuCn  +  Z.vijUzjCO 


(16) 


where  U  is  the  internal  energy  of  the  phonon  system.  Although 
Eq.  (16)  gives  the  correct  temperature  dependent  trend  assuming 
that  all  phonon  frequencies  in  each  material  available  at  any  given 
temperature  participate  equally  to  energy  transmission,  it  does  not 
differentiate  between  inelastic  and  elastic  contributions  nor  does  it 
take  into  account  the  relative  contributions  from  each  phonon  fre¬ 
quency  (i.e.,  it  assumes  that  all  phonon  frequencies  have  an  equal 
probability  of  participating  in  phonon  scattering  events).  How¬ 
ever,  there  is  some  intrinsic  frequency  dependence  on  the  trans¬ 
mission  coefficient,  which  is  contained  in  the  temperature  depen¬ 
dence  of  the  experimental  data  and  not  contained  in  the 
assumptions  of  Eq.  (16). 

In  order  to  determine  the  relative  contributions  from  elastic  and 
inelastic  scattering,  the  coefficients  A  and  6(7)  must  be  examined 
in  greater  detail,  which  can  be  accomplished  by  examining  hBD 
data  clearly  dominated  by  inelastic  scattering.  In  this  analysis, 
only  temperatures  in  the  classical  limit  of  the  metal  films  will  be 
considered  (i.e.,  temperatures  above  the  softer  material’s  Debye 
temperature).  In  this  limit  the  data  show  the  linear  trend  associ¬ 
ated  with  inelastic  scattering;  however,  a  fraction  of  this  hBD  is 
expected  to  be  associated  with  elastic  scattering,  which  is  ex¬ 
pected  to  be  constant. 


To  determine  hBD,  the  linear  fit  to  the  data  was  extended  to  the 
y-axis  to  find  the  y-intercept,  which  is  taken  as  the  contribution  of 
elastic  scattering  in  the  classical  limit.  The  linear  fit  of  the  Pb/ 
diamond  data  compared  to  /7g^L,  hBD,  and  h™ M  is  shown  in  Fig. 
5.  The  y-intercept  is  hBD  =  3.26  MW  m-2  K-1  where 
=  6.17  MW  m-2  K-1.  As  expected,  /2^,L  is  greater  than  hBD  deter¬ 
mined  from  experimental  data,  since  not  all  substrate  phonons  will 
be  participating  in  conductance.  Comparing  these  two  values 
yields  A  =  0.53.  In  this  limit,  /i™m=3.32  MW  m-2  K-1,  which  is 
in  excellent  agreement  with  hBD,  indicating  that  the  DMM  pre¬ 
dicts  the  elastic  contributions  to  hBD  at  the  heavily  mismatched 
Pb/diamond  interface  well.  Table  2  compares  values  of  heBD  (de¬ 
termined  from  the  y-intercept),  and  the  coefficient  A 

for  six  samples  in  the  classical  limit.  Also  listed  in  Table  2  is  the 
ratio  of  Debye  temperatures  of  the  materials  comprising  side  1 
and  side  2  for  each  interface.  This  ratio  quantifies  the  degree  of 
acoustic  mismatch  of  each  interface — the  smaller  the  ratio,  the 
greater  the  mismatch  between  phonon  spectra.  Note  that  the 
agreement  between  heBD  and  is  much  better  for  the  heavily 

mismatched  samples  than  for  the  better  matched  Pt  samples.  This 
is  expected  since  the  transmission  coefficient  calculations  become 
less  accurate  as  the  materials  become  more  similar  [23],  This  in¬ 
dicates  that  the  DMM  is  a  good  model  for  predicting  the  elastic 
scattering  contribution  to  thermal  boundary  conductance,  and  the 
transmission  coefficient  calculated  with  Eq.  (5)  becomes  more 
accurate  as  the  adjacent  phonon  spectra  become  more  dissimilar. 

Now  that  the  elastic  portion  of  thermal  boundary  conductance 
has  been  determined,  the  inelastic  contribution  can  be  estimated 
from  the  slope  of  the  experimental  data.  Rearranging  Eq.  (14),  the 
coefficient  representing  the  transmission  of  substrate  phonons 
through  inelastic  scattering  can  be  expressed  as 


B(T)  = 


hBD(T)-AhpB g- 


h 


IPRL 

BD 


(17) 


where  hBD{T)  is  the  thermal  boundary  conductance  as  a  function 
of  temperature  determined  from  the  linear  fit  to  the  experimental 
data,  and  A  was  determined  above.  Note  that  from  the  discussion 
of  Eq.  (15),  B(T)  =  a'2e\T).  The  coefficient  B(T)  as  a  function  of 
temperature  can  be  calculated  only  over  the  temperature  range  of 
the  experimental  data.  However,  higher  temperature  values  were 
determined  through  a  nonlinear  regression  extrapolation  to  deter¬ 
mine  6(7)  at  all  temperatures  T>  eD.  Due  to  the  nature  of  this 
high  temperature  extrapolation,  above  a  given  temperature  the 
slope  of  the  IPRL  decreases  to  a  value  that  is  less  than  the  slope  of 
hBD(T),  which  causes  an  unphysical  slight  linear  increase  in  6(7) 
resulting  in  a  local  minimum.  To  correct  for  this  unphysical  in¬ 
crease  in  Eq.  (17),  at  temperatures  above  the  local  minimum  tem¬ 
perature,  6(7)  is  fixed  to  a  constant  value  of  the  local  minimum. 
This  gives  a  much  more  physical  trend  to  6(7)  and  avoids  any 
error  due  to  the  nonlinear  extrapolation  routine.  Figure  6  com¬ 
pares  the  temperature  dependent  transmission  coefficient  that 


Table  2  High  temperature  limits  of  thermal  boundary  conductance  and  the  elastic  contribution 
to  hBD,  the  DMM,  the  PRL,  and  the  predicted  hBD  based  on  experimental  trends  [hBD(7)]  for  six 
different  interfaces.  All  units  are  in  MWm_2K_1,  except  A,  B,  and  0D1/0D, 2>  which  are  unitless 
ratios. 


Interface 

0D,l/ @D,2 

/iel 

nBD 

i  DMM 
nBD 

hPRL 

nBD 

A 

B 

/!bd(7=°°) 

Pb/diamond 

0.047 

3.26 

3.31 

6.23 

0.523 

0.00436 

58.7 

Pb/H/diamond 

0.047 

1.92 

3.31 

6.23 

0.310 

0.00329 

41.9 

Bi/H/diamond 

0.053 

2.25 

1.32 

1.34 

1.67 

0.00158 

22.3 

Au/diamond 

0.074 

16.6 

14.4 

31.2 

0.532 

0.00626 

96.1 

Pt/Al203 

0.21 

51.8 

130 

190 

0.261 

0.0170 

210 

Pt/AIN 
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Fig.  6  Comparison  of  a4nel(7)  determined  by  1-a!,nel(7)  where 
a!,nel(7)  is  given  by  the  model  derived  by  Dames  and  Chen  [9] 
and  presented  in  Eq.  (16);  and  6(7),  which  is  derived  in  this 
work  and  presented  in  Eq.  (17).  6(7)  decreases  to  a  much  lower 
constant  value  much  more  quickly  than  afel(T)  calculated  with 
Eq.  (16).  The  assumption  in  Eq.  (16)  is  that  all  phonons  of  all 
frequencies  have  equal  probability  in  transmission,  where  Eq. 
(17)  is  based  on  experimental  data  and  the  IPRL.  The  smaller 
phonon  transmission  of  side  2  phonons  at  higher  temperatures 
could  be  due  to  the  probability  of  side  2  phonons  breaking 
down  into  lower  frequency  side  1  phonons  decreasing  as  the 
side  2  phonon  frequency  increases. 


takes  into  account  inelastic  scattering  for  side  2  derived  by  Dames 
and  Chen  [9],  a'i,nel(T)=  1  -c^el(T),  where  a'l1nel(7')  is  calculated 
with  Eq.  (16),  to  the  inelastic  transmission  coefficient  determined 
from  Eq.  (17),  B(T),  for  a  Pb/diamond  interface.  B(T)  decreases  to 
a  much  lower  constant  value  much  more  quickly  than 
calculated  with  Eq.  (16).  As  previously  mentioned,  Eq.  (16)  as¬ 
sumes  that  all  phonon  frequencies  participate  equally  in  transmis¬ 
sion  and  hBD.  B(T),  however,  takes  into  account  the  temperature 
dependent  phonon  scattering  trends  in  inelastic  processes  by  uti¬ 
lizing  the  trends  encapsulated  in  the  experimental  data.  The  values 
and  trends  in  B(T)  show  that  as  temperature  increases  and  more 
high  frequency  phonons  become  available  in  side  2  that  are  not 
available  in  side  1  (i.e.,  the  phonon  frequencies  in  side  2  that  are 
above  the  side  1  cutoff  frequency),  the  transmission  of  energy 
from  side  2  to  side  1  is  much  less  than  if  assuming  all  phonon 
frequencies  participate  equally  in  energy  transmission  as  assumed 
in  Eq.  (16).  This  suggests  that  phonons  in  side  2  with  frequencies 
much  higher  than  the  side  1  cutoff  are  less  likely  to  participate  in 
inelastic  processes  than  phonons  with  frequencies  that  are  only 
slightly  higher  than  the  side  1  cutoff.  A  similar  trend  was  inferred 
by  Hopkins  et  al.  [17],  and  was  attributed  to  the  probability  a 
three-phonon  process  occurring,  which  could  occur  for  coupling 
between  a  lower  frequency  side  2  phonon  and  two  side  1  phonons, 
being  higher  than  the  probability  of  ;?-phonon  processes  occurring 
(where  n> 3),  which  would  have  to  occur  for  a  higher  frequency 
side  2  phonon  to  couple  with  n  - 1  side  1  phonons. 

With  this,  the  contributions  of  elastic  and  inelastic  scattering  to 
hBD  were  calculated  based  on  experimental  data  and  extended  to 
high  temperatures  to  determine  hBD  in  the  classical  limit  that  takes 
into  account  inelastic  scattering.  The  inelastic  and  elastic  contri¬ 
butions  to  hBD  for  Pb/diamond  are  shown  in  Fig.  7,  along  with 
1iBd{T)  calculated  from  Eq.  (14).  The  high  temperature  trends  fol¬ 
low  what  is  expected  when  temperatures  are  driven  higher  than 
both  materials’  Debye  temperatures.  This  same  agreement  be- 
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Fig.  7  Inelastic  and  elastic  scattering  contributions  to  hBD  of 
Pb/diamond  in  the  classical  limit.  The  total  thermal  boundary 
conductance,  /jbd(7),  shows  excellent  agreement  to  the  data  at 
low  temperatures  in  the  classical  limit,  and  approaches  a  con¬ 
stant  value  at  higher  temperatures  that  are  above  the  Debye 
temperature  of  diamond  (or  at  temperatures  that  are  greater 
than  the  Debye  temperature  of  both  materials). 


tween  Eq.  (14)  and  the  experimental  data  and  high  temperature 
trends  are  shown  between  Eq.  (14)  and  the  Bi/diamond,  Au/ 
diamond,  Pt/Al203,  and  Pt/AIN  data. 

To  demonstrate  the  relative  magnitude  of  inelastic  scattering  on 
hBD,  the  ratio  h'BplheBD  is  plotted  versus  temperature  in  the  clas¬ 
sical  regime  in  Fig.  8  for  the  six  interfaces  studied.  As  expected, 
the  inelastic  contribution  compared  to  the  elastic  contribution  to 
hBD  in  Pb,  Bi,  and  Au  on  diamond  are  more  temperature  depen¬ 
dent  than  the  relative  contribution  in  the  Pt/Al203  and  Pt/AIN 
samples  due  to  the  higher  Debye  temperature  of  diamond  than 
A1203  and  AIN.  Also,  the  contribution  of  inelastic  scattering  to 
overall  hBD  is  greater  than  that  of  elastic  scattering,  and  the  rela¬ 
tive  contribution  of  inelastic  scattering  increases  with  interface 


Fig.  8  Relative  magnitude  of  inelastic  scattering  on  hBD.  This 
ratio  compares  the  inelastic  contribution  to  hBD  to  the  elastic 
contribution  as  predicted  via  Eq.  (14).  The  role  of  inelastic  pho¬ 
non  scattering  increases  as  the  acoustic  mismatch  of  the  film 
and  substrate  becomes  greater.  The  range  in  which  hBD  should 
increase  linearly  with  temperature  due  to  inelastic  scattering 
also  increases  with  acoustic  mismatch. 
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acoustic  mismatch  (i.e.,  as  the  film/substrate  Debye  temperature 
ratio  become  smaller).  This  also  leads  to  the  inelastic  contribution 
becoming  independent  of  temperature  at  higher  temperatures  in 
the  diamond  samples  than  in  the  Pt  on  A1203  and  AIN  samples.  In 
fact,  the  high  temperature  limit  of  the  Pb,  Bi,  and  Au  on  diamond 
samples  is  not  reached  until  the  interface  temperature  is  driven 
above  the  melting  temperature  of  the  metal,  meaning  that  in  nano¬ 
structures  with  Pb.  Bi,  and  Au  films  adjacent  to  diamond  struc¬ 
tures,  hBD  will  continually  increase  with  temperature  until  melt¬ 
ing.  Also,  in  the  high  temperature  limit,  as  the  mismatch  between 
the  materials  adjacent  to  the  interface  grows,  the  contribution  of 
inelastic  phonon  scattering  to  thermal  boundary  conductance  in¬ 
creases. 

Conclusions 

The  accuracy  of  hBD  predictions  made  using  DMM  and  PRL 
varies  depending  on  the  types  of  materials  comprising  the  inter¬ 
face.  Due  to  inelastic  scattering  these  models  may  not  be  valid  at 
temperatures  characteristic  of  many  modern  nanodevices  and 
structures  since  these  models  do  not  take  into  account  inelastic 
phonon  scattering.  Although  the  JFDMM  provides  a  simple  ap¬ 
proximation  of  the  effects  of  joint  vibrational  modes  on  hBD  using 
a  Debye  approximation,  it  does  not  specially  examine  the  contri¬ 
butions  due  to  elastic  and  inelastic  scattering.  These  contributions 
can  be  extracted  from  the  trends  and  values  in  the  experimental 
data.  To  estimate  the  role  of  elastic  scattering  in  the  classical  limit, 
the  v-intercept  of  the  linear  trend  in  the  experimental  data  of  hBD 
as  a  function  of  temperature  is  compared  to  the  PRL  in  the  clas¬ 
sical  limit.  The  resulting  contribution  of  elastic  scattering  in  this 
limit  agrees  well  with  the  value  of  hBD  predicted  by  the  DMM. 
This  agreement  increases  with  sample  mismatch,  which  validates 
the  DMM  as  a  model  to  predict  the  contribution  of  elastic  scatter¬ 
ing  in  heavily  mismatched  samples.  To  estimate  the  role  of  inelas¬ 
tic  scattering  in  the  classical  limit,  the  slope  of  the  experimental 
data  of  hBD  as  a  function  of  temperature  is  compared  to  a  new 
model,  the  IPRL.  The  IPRL  predicts  the  maximum  thermal  bound¬ 
ary  conductance  assuming  all  substrate  phonons  of  all  frequencies 
are  transmitted  into  the  film.  The  IPRL  predictions  are  compared 
to  the  experimental  data  to  effectively  determine  what  percentage 
of  substrate  phonons  are  participating  in  inelastic  scattering  be¬ 
yond  those  participating  in  elastic  scattering.  The  role  of  inelastic 
scattering  to  overall  hBD  is  shown  to  outweigh  the  role  of  elastic 
scattering  at  the  interfaces  of  interest,  with  an  increasing  role  of 
inelastic  scattering  with  interfacial  acoustic  mismatch.  The  predic¬ 
tions  of  inelastic  scattering  show  a  linear  trend  at  lower  tempera¬ 
tures  in  the  classical  limit  and  predict  a  constant  hBD  at  higher 
temperatures,  giving  an  upper  limit  to  thermal  boundary  conduc¬ 
tance  that  takes  into  account  both  elastic  and  inelastic  scattering. 
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Nomenclature 

A  —  coefficient  relating  to  elastic  scattering  in  Eq. 

(14) 

B  =  coefficient  relating  to  inelastic  scattering  in  Eq. 

(14) 


D  =  phonon  density  of  states  per  unit  volume, 
s  m-3 

fi  =  Planck’s  constant  divided  by  2tt,  J  s 
h BD  =  thermal  boundary  conductance,  W  m-2  K-1 
M  =  molecular  weight,  g  mol-1 
N  =  phonon  number  density,  m-3 
Na  =  Avogadro's  number,  mol-1 
n  =  Bose-Einstein  distribution  function 
q  =  heat  flux,  W  m-2 
T  =  temperature,  K 
U  =  internal  energy,  J  m-3 
v  =  phonon  group  velocity,  m  s_1 


Greek  Symbols 

a  = 

$  = 
w  = 

£  = 

Subscripts 

1  = 

2  = 

j  = 
L  = 
mod  = 
T  = 

Superscripts 

c  = 

DMM  = 
el  = 
inel  = 
IPRL  = 

PRL  = 
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elastic 

inelastic 

calculated  with  the  inelastic  phonon  radiation 
limit 
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